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We study x-ray absorption spectra of azobenzene-functionalized self-assembled monolayers 
(SAMs), investigating excitations from the nitrogen K edge. Azobenzene with H-termination and 
functionalized with CFg groups is considered. The Bethe-Salpeter equation is employed to compute 
the spectra, including excitonic effects, and to determine the character of the near-edge resonances. 
Our results indicate that core-edge excitations are intense and strongly bound: Their binding ener¬ 
gies range from about 6 to 4 eV, going from isolated molecules to densely-packed SAMs. Electron- 
hole correlation rules these excitations, while the exchange interaction plays a negligible role. 

PACS numbers: 71.35.-y, 73.20.Mf, 78.70.Dm 


I. INTRODUCTION 

Azobenzene-functionalized self-assembled monolayers 
(SAMs) on metal surfaces represent a viable and efficient 
way to obtain ordered architectures of photo-switching 
moleculesHowever, it has been observed that in such 
closely packed systems photo-isomerization can be dras¬ 
tically hindered by steric effects,!^®^ and even by exci¬ 
tonic coupling between the chromophores.^ In order to 
overcome these limitations and obtain SAMs with effi¬ 
cient switching rates, a number of strategies have been 
developed, such as modifying the mor pholo gy of the 
substrate,^ introducing organic spacersp^ and func¬ 
tionalizing azobenzene with end groups. 

To tune these complex systems in view of optimized 
performance, a deep knowledge of their chemical com¬ 
position and structure-property relationship is required, 
x-ray absorption spectroscopy (XAS) represents a pow¬ 
erful technique for this purpose, and a synergistic inter¬ 
play with theory can provide an insightful interpretation 
of the experimental data. First-principles methods rep¬ 
resent the most suitable tool. Density-functional theory 
(DFT), both in the core-hole approximatiorP^ and the A- 
self-consistent-field (ASCF) approach,!^ is routinely ap¬ 
plied to simulate XAS in a wide range of materials, from 
gas-phase molecules to solid-state systemsRecently, 
also time-dependent DFT has become popular to com¬ 
pute core-level excitations in molecular compounds.®^ 
While these approaches can provide qualitative agree¬ 
ment with experiments, explicit many-body treatment 
has turned out superior to such approaches. In small 
molecules, coupled-cluster methods have been success¬ 
fully applied t o com pute XAS from the carbon and ni¬ 
trogen K-edge.l^^J^ For solid-state materials, many-body 
perturbation theory (MBPT) represents the st ate-of -the- 
art formalism to describe neutral excitations) ^^ * ^^ * The 
electron-hole {e-h) interaction, effectively described by 
the Bethe-Salpeter equation (BSE), pla ys a c rucial role 
not only in conventional semiconductorsp^^ but also in 
organic crystalP^®^ and even in isolated molecules .^2 A 
number of studies on core-level excitations, from differ- 
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EIG. 1: (Color online) a) Sketch of azobenzene-functionalized 
SAM of alkanethiols on a gold surface, b) isolated azobenzene 
molecule, c) diluted (d-SAM) and d) packed SAM (p-SAM) in 
their unit cells, including one and two inequivalent molecules, 
respectively. 


ent edges and in several materials,E^®^has demonstrated 
that BSE can accurately reproduce XAS. 

In this paper, we present an ab initio study of x-ray 
absorption spectra of azobenzene-functionalized SAMs. 
We consider excitations from the nitrogen (N) K edge, 
i.e., involving transitions from Is electrons to the conduc¬ 
tion bands. In this manner, we obtain exciton binding 
energies and determine the character of the core-level ex¬ 
citations. Going from the isolated molecule to a closely- 
packed SAM, we analyze the XAS at increasing density of 
azobenzene molecules, and we compare our results with 
experimental data. In order to understand whether and 
how functional groups affect the nature of the excitons 
and their binding energy, we consider molecules that are 












2 


either H-terminated (H-az) or functionalized with triflu- 
oromethyl (CFg-az). 


II. SYSTEM AND METHODS 


A. Azobenzene SAMs 


A sketch of an azobenzene-functionalized SAM of alka- 
nethiols on gold is presented in Fig. The chro- 

mophores are covalently bonded to the alkyl chains, 
which are attached to the gold substrate through a thiol 
group. As suggested by s canning tunneling microscopy 
(STM) measurement s,^i^ the SAM has an orthorhombic 
supercell, with lattice vectors a=6.05 A and 6=7.80 A, 
hosting two inequivalent azobenzene molecules. In our 
calculations, we neglect the alkyl chains and the gold 
surface, since they are expected not to play a role in the 
XAS from the N K edge. Therefore, we consider only the 
azobenzene molecule, with a methoxy group added to one 
end (see Fig. in order to reproduce the chemical en¬ 
vironment of the covalent bond to the alkyl chain. The 
reciprocal distance and orientation of the molec ules in the 
unit cell is set according to the STM data.^^^^ Although 
the first experiments on these systems predicted a her¬ 
ringbone structure of the chromophores in the SAMs,l^ 
a consensus about the orientation of azobenzenes is still 
missing. We consider the two inequivalent molecules in 
the unit cell being oriented parallel to each other, since 
we expect deep core levels to be hardly affected by the 
reciprocal orientation of the molecules. In this configu¬ 
ration (see Fig. Et) the azobenzenes are separated by 
about 2 A in the lateral direction, and by ^ 3.8 A in the 
direction perpendicular to the plane of the phenyl rings. 
We incorporate ^ 14 A of vacuum in the vertical direc¬ 
tion, to effectively simulate a two-dimensional system. 

In order to understand the effects of packing in the 
XAS, we consider an additional structure, including only 
one molecule in the same unit cell. We refer to this sys¬ 
tem, shown in Fig. [^, as diluted SAM (d-SAM), to 
distinguish it from the packed SAM (p-SAM, Fig. Et)- 
For comparison, we investigate an isolated azobenzene 
molecule in an orthorhombic supercell, with ^ 6 A of 
vacuum in each direction. We also consider SAMs of 
CFg-functionalized azobenzene, which have been recently 
synthesized.fiS CF 3 and other functional groups are used 
in experiments as markers, to identify the orientation of 
the molecules with respect to the surface,!^ and/or to 
tune the switching properties of the SAMs by decreasing 
the steric hindrance due to inter molecular interactions.^ 
Also for CFg-az, we investigate p- and d-SAMs, as well 
as an isolated molecule for comparison. We adopt the 
same structures shown for H-az in Fig. 


B. Theoretical Background 

x-ray absorption spectra are computed from first prin¬ 
ciples by solving the BSE, which is an effective equa¬ 
tion of motion for the electron-hole two-particle Green’s 
fnnctionj ^^ l ^^ l By considering only transitions from core 
(c) to unoccupied (u) states, the BSE in matrix form 
reads: 


E ttBSE aX _ z?A jA 

^cnk,c'u'k'^c'w'k' “ ^ ^cuk* 

c'u'k' 


( 1 ) 


In case of N A edge, the N Is is the only initial state. 
The BSE Hamiltonian in Eq. can be written as the 
sum of three terms: 


^BSE ^ fydiag ^ . (2) 


The diagonal term accounts for single-particle 

transitions. Including only this term corresponds to the 
independent-particle approximation (IPA). The exchange 
(H^) and direct terms in Eq. [^incorporate the 

repulsive bare and the attractive screened Coulomb in¬ 
teraction, respectively. The coefficients 7 ^, and 7 c in Eq. 
[^ enable to select the spin-singlet ( 7 ^, = 7 c = 1 ) and 
spin-triplet ( 7 ^, = 0, 7 c = 1) channels. In the latter 
case, the exchange interaction is not present. In Eq. 
the eigenvalues represent excitation energies. Exciton 
binding energies (E 5 ) are defined, for each excitation, as 
the difference between excitation energies computed 
from IPA and BSE. The eigenvectors carry informa¬ 
tion about the character and composition of the excitons. 
Through the transition coefficients 


Ia 


E 

cwk 


^cwk 


(ck|p|nk) 

^uk ^ck 


(3) 


enter the expression of the imaginary part of the 
macroscopic dielectric function (^m)- 

ImsM = - Ex). (4) 

A 


C. Computational Details 

All calculations are performed with the exciting 
code,l^ a computer package implementing DET and 
MBPT.I^ exciting is based on the all-electron full- 
potential augmented planewave method, which ensures 
an explicit and accurate description of core electrons. 
The calculation of XAS via the solution of the BSE 
in an all-electron framework has been successfully ap¬ 
plied to different abs orption edges in a number of 
bulk materials,S^®^EHl including, very recently, small 
molecules . 1 ^ 

The Kohn-Sham (KS) electronic structure, used 
here as starting point for the BSE, is computed 
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within the local-density approximation (Perdew-Wang 
functional)!^ Quasiparticle energies are approximated 
by KS single-particle energies, and a scissors operator 
is applied to match the experimental absorption onset 
for the p-SAM, according to the available data for H- 
aJEl and CFg-az.f^ The same correction is applied also 
to the respective isolated molecule and d-SAM, since 
we do not have experimental data available for these 
systems!^ A k-point mesh of 6x4x1 (3x2x1) is used 
to sample the Brilloiun zone for the p-SAM (d-SAM), in 
both ground-state and BSE calculations. For the basis 
functions, a planewave cutoff Gmax=^ bohr is applied 
to the molecules; for the SAMs it is reduced to 4.625 
bohr. Mufhn-tin spheres of radii i?MT=0.8 bohr are con¬ 
sidered for hydrogen, bohr for nitrogen and 

fluorine, and i?MT=1-2 bohr for carbon and oxygen. The 
atomic positions of each structure are optimized by min¬ 
imizing the Hellmann-Feynman forces within a threshold 
of 0.025 eV/A. For the calculation of the e-h interaction 
term in the BSE (Eq. [^, the screening is evaluated 
within the random-phase approximation, including the 
N Is core states, all valence states and 100 conduction 
bands. Local-field effects are included, with at least 400 
|G + q| vectors for the SAMs and about 2000 for the 
molecules. These parameters ensure convergence of the 
XAS within 0.25 eV. 


III. RESULTS AND DISCUSSION 

In Fig. we present the XAS computed for the 
azobenzene molecule and for the p- and d-SAMs. In ad¬ 
dition to the result obtained from the solution of the 
BSE, the corresponding IPA spectrum is shown in each 
panel. To better guide the reader, we start our analy¬ 
sis from the IPA results. These are presented only for 
comparison with the BSE spectra, to better highlight ex- 
citonic effects. Three peaks, labeled as A, B and C, can 
be identified approximateiy at the same energy in each 
structure, independentiy of the packing density. Since 
each moiecuie has two N atoms, with the Is ieveis sep¬ 
arated by 61 meV in the KS spectrum, this muitipiicity 
is reflected also in the XAS. Hence, in the spectrum of 
the isolated molecule {top panel) each peak is formed 
by two transitions. The first and most intense peak. A, 
corresponds to a transition to the LUMO ievei, which 
presents tt* character, with strong iocaiization on the 
azo group (see Fig. [^. Aiso the other peaks, B and 
C, are given by tt* resonances, invoiving excitations to 
LUMO-l-3 and LUMO+7, respectiveiy. By inspecting 
the orbitals, shown in Fig. we observe a direct cor¬ 
respondence between the amount of charge localized on 
the azo group and the strength of the resonance. Peaks 
A and C have large intensity, while B is rather weak. Be¬ 
tween A and B additional transitions to the unoccupied 
levels LUMO+1 and LUMO+2 take place. Considering 
that these orbitals are localized predominantly on the 
phenyl rings (Fig. |^, the corresponding transitions are 
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FIG. 2: (Color online) XAS of isolated azobenzene molecule 
(top), d-SAM (middle), and p-SAM (bottom). In the BSE 
spectra (solid line), the peaks are labelled as A’, B’ and C’. 
Independent-particle spectra (IPA - solid area) are shown 
for comparison (peaks A, B and C). The imaginary part of 
the macroscopic dielectric function is averaged over the three 
Cartesian components. The experimental curve for the p- 
SAM is taken from Ref. [151 A Lorentzian broadening of 0.1 
eV is applied to the theoretical spectra. H-az is shown on top. 


extremely weak. 

The independent-particle picture for the SAMs 
presents analogous features observed for the isolated 
molecule. Interestingly, not only the relative intensity 
of the three peaks is the same, but also their energy is 
independent of the packing density. In the SAMs a larger 
number of transitions contribute to the peaks, compared 
to the single molecule. This is especially true for the p- 
SAM, where the peaks A and C experience a broadening 
of about 1 eV and 0.5 eV, respectively (see Table In 
the spectrum of the d-SAM, 8 single-particle transitions 
contribute to the lowest-energy peak A, which is com¬ 
prised within less than 0.1 eV. Also the nature of the 
transitions is the same in the single molecule and in the 
SAMs. 

The scenario changes significantly when e-h interaction 
is taken into account. By inspecting the BSE spectrum 
of the single molecule (Fig. top panel), we observe 
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Method 

BSE 

IPA 

Peak 

A’ 

B’ 

C’ 

A 

B 

C 

Molecule 

396.3 [ 6 . 0 ] 

400.8 [ 3 . 7 ] 

401.2 [ 4 . 2 ] 

402.3 

404.5 

405.4 

d-SAM 

396.8 [ 5 . 7 ] 

401.3 [ 3 . 3 ] 

401.7 [3.6] 

402.5 

404.6 - 404.8 

405.3 - 405.6 

p-SAM 

398.0 [4.1] 

402.0 [1.9] 

402.4 [ 2 . 9 ] 

402.1 - 402.9 

403.9 - 405.0 

405.3 - 405.9 


TABLE I: Excitation energies of the main peaks of the spectra presented in Fig. for molecule and SAMs. A scissors operator 
of 26.8 eV is applied to the underlying KS electronic structures to match the experimental absorption onset (Ref. IThll . Exciton 
binding energies, corresponding to the difference between IPA and BSE excitation energies (i.e., Ea — Ea', Eb — Eb', and 
Ec — Eci), are highlighted in bold. All quantities are expressed in eV. 
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FIG. 3: (Color online) Lowest unoccupied molecular orbitals 
of the azobenzene molecule. An isovalue of 0.002 is used for 
the isolsurfaces. 


again three peaks, labeled as A’, B’ and C’, whose ener¬ 
gies are significantly red-shifted compared to their IPA 
counterparts (see Table The exciton corresponding to 
the peak A’ has a sizable binding energy of 6.0 eV. Also 
B’ and C’ are remarkably red-shifted compared to their 
IPA counterparts, by 3.7 eV and 4.2 eV, respectively. In¬ 
terestingly, the e-h interaction is stronger in C’ than in 
B’. It is worth noting that the relative energy difference 
between A’ and B’ (4.6 eV) is more than twice than the 
one between A and B (2.2 eV). Also the energy difference 
between A’ and C’ (4.9 eV) exceeds by about 50% the 
one between A and C (3.1 eV). In contrast, B’ and C’ are 
separated only by 0.4 eV, which is half of the difference 
between B and C (0.9 eV). Also the oscillator strengths 
are considerably redistributed upon inclusion of the e- 
h interaction. A’ is almost twice more intense than A, 


with a consequent decrease of relative spectral weight of 
B’ and C’. From the analysis of the exciton composition, 
we observe that A’ corresponds to a pure transition to 
the LUMO, as within the IPA. Conversely, B’ presents 
a rather mixed character, with a dominant contribution 
(^70%) from the transition to the LUMO+3 state. Also 
C’ is composed by transitions to LUMO+7 (^ 70%) and 
LUMO (-12%). 

The main features observed in the XAS of the single 
molecule can be recognized also for the SAMs (Fig. 
middle and bottom panels). The lowest-energy peak A’ 
dominates the spectrum in terms of spectral weight, and 
its energy is red-shifted with respect to A by 5.7 and 4.1 
eV in the d- and p-SAM, respectively (see also Table |I|. 
In the case of the d-SAM, the peaks B’ and C’ are signifi¬ 
cantly red-shifted with respect to their IPA counterparts 
B and C, by 3.3 and 3.6 eV, respectively. Like for the 
molecule, also in this case the binding energy of C’ is 
larger than the one of B’. In the spectrum of the p-SAM, 
— 3 eV for C’, and < 2 eY for B’ (see Table |I|. 
These values of binding energies, especially for the core¬ 
edge exciton A’, are remarkably larg^^ with respect to 
optical excitations in molecular crystals, where typically 
Eb < 1 eV.I^ They are also large compared to inorganic 
materials: For example, near-edge resonances from the Li 
or Be K edge in binary crystals present binding en ergies 
ranging from a few hundreds of meV up to — 2 

The BSE spectra in Fig. are almost rigidly blue- 
shifted of about 2 eV, when going from the isolated 
molecule to the p-SAM. Considering that the IPA absorp¬ 
tion onset, given by the position of A, is approximately 
the same in the three systems (see also Table |^, the de¬ 
creasing binding energy of A’ upon increasing azobenzene 
packing density is driven by many-body interactions. We 
assign this effect to a combination of enhanced screening 
and wave-function delocaliza tion, as observed in optical 
excitations of organic crystals.^^^®^ Moreover, as a con¬ 
sequence of dipole-dipole coupling, the plethora of single¬ 
particle transitions forming the lowest-energy peak in the 
SAMs, combine into two excitons, corresponding to the 
peak A’P 

The predominant role of the e-h correlation in such 
core-edge excitons is further confirmed by an additional 
analysis of the BSE results. When computing triplet ex¬ 
citation energies, we diagonalize a BSE Hamiltonian with 
lx = 0 (Eq. §, since triplet states do not experience 
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exchange interaction. By comparing singlet and triplet 
excitation energies, we notice differences of the order of 
20 - 50 meV for all the considered excitons. This result 
is independent of the packing density of the molecules 
and reveals that the screened Coulomb interaction is the 
driving mechanism of the XAS. Conversely, local-field ef¬ 
fects (LFE), ruled by the exchange term ^ are irrele¬ 
vant here, due to the localized character of the excitation. 
There is in fact a very small overlap between the h, lo¬ 
calized in the N Is state, and the e in the conduction 
region. This represents a rather different scenario than 
optical absorption, where LFE are significantly larger. In 
particular, organic materials are typically characterized 
by singlet-triplet splitting of t he same o rder of magni¬ 
tude of singlet binding energiesConcerning the 
character of the excitons in the SAMs, the picture is not 
significantly different from that of isolated azobenzene. 
The exciton A’ corresponds to transitions from N Is to 
the lowest unoccupied band. In the p-SAM, which has 
two inequivalent molecules in the unit cell, this band is 
split in two. Transitions to both subbands mix up to 
form the exciton. The peak B’ is given by a manifold of 
weak excitations, with a rather mixed character, in both 
d- and p-SAM. They include transitions to unoccupied 
states, which are the counterparts of LUMO+3 in the 
single molecule (see Fig. |^. The peak C’ is composed by 
only 2 and 4 excitons in the d- and p-SAM, respectively. 
These excitons have a remarkably mixed character, and 
involve transitions to higher unoccupied bands. 

The quality of our BSE results is confirmed by the 
comparison with available experimental data for the p- 
SAMBy inspecting Fig. i we notice that the ex¬ 
perimental spectrum is characterized by a strong peak, 
corresponding to A’. It is attributed to a tt* resonance, 
due to the Is LUMO transition?!^ this is in perfect 
agreement with our finding. About 4 eV above, a shoul¬ 
der appears in correspondence of B’ and C’. Due to the 
limited resolution of the experimental spectrum, it is not 
possible to identify these two excitons separately. How¬ 
ever, also in the experiment, they are assigned to transi¬ 
tions to LUMO-Fn states, having again tt* character,!!^ as 
confirmed by our results. It is worth noting that besides 
the peaks A’, B’ and C’, the experimental spectrum has 
non-zero intensity around 401 eV. This feature is present 
also in our BSE spectrum with extremely small oscilla¬ 
tor strength.l^ This weak peak is given by two double 
degerate excitations, targeting the two lowest unoccu¬ 
pied bands (LUMO and LUMO-fl). Furthermore, in the 
high-energy part of the experimental spectrum, the oscil¬ 
lator strength increases towards ^ 408 eV . w here a bump 
due to the Is ^ a* transition is observed.^ This feature 
is not reproduced by theory, since cr* states, which lie 
very high in energy in the KS spectrum, are not consid¬ 
ered in the solution of the BSE. 

Finally, we investigate the x-ray absorption spectra 
of azobenzene molecules and SAMs terminated with a 
trifluoro-methyl group (CFg-az, see Fig. [^. Since no 
additional N atoms are introduced in the system, we ex- 
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FIG. 4: (Color online) XAS of isolated CFg-functionalized 
azobenzene molecule (top), d-SAM (middle), and p-SAM 
(bottom). In the BSE spectra (solid line), the peaks are la¬ 
belled as A’, B’ and C’. Independent-particle spectra (IPA - 
solid area) are shown for comparison (peaks A, B and C). The 
imaginary part of the macroscopic dielectric function is aver¬ 
aged over the three Cartesian components. The experimental 
curve for the p-SAM is taken from Ref. M A Lorentzian 
broadening of 0.1 eV is applied to the theoretical spectra. 
CFg-az is shown on top. 


pect to observe the same features as in the XAS of the 
H-terminated counterparts. With this analysis we aim at 
understanding the influence of functionalization on bind¬ 
ing energies and exciton character. We again consider an 
isolated CFg-functionalized azobenzene molecule, as well 
as d- and p-SAMs, in order to inspect the role of packing 
density. The calculated XAS are shown in Fig. and 
the (binding) energies of the bright excitons are summa¬ 
rized in Table [llj The spectra appear strikingly similar 
to those presented in Fig. for H-terminated azoben¬ 
zene, and so the main features analyzed previously. The 
XAS are considerably blue-shifted from the molecule to 
the p-SAM, i.e., upon increasing intermolecular interac¬ 
tions. The intense peaks correspond to transitions to 
unoccupied states, having the same tt* character as in 
the H-az systems. By inspecting carefully Tables |I] and 
|H| we notice that exciton binding energies slightly in- 
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Method 

BSE 

IPA 

Peak 

A’ 

B’ 

C’ 

A 

B 

C 

Molecule 

396.6 [6.3] 

401.1 [3.9] 

401.6 [4.4] 

402.9 

405.0 

406.0 

d-SAM 

397.5 [5.5] 

401.7 [3.5] 

402.2 [3.7] 

403.0 

405.2 

405.9 - 406.1 

p-SAM 

398.3 [4.5] 

402.5 [2.1] 

403.0 [2.8] 

402.8 - 403.4 

404.6 - 405.7 

405.8 - 406.4 


TABLE II: Excitation energies of the main peaks of the spectra presented in Eig. for CEg-az molecule and SAMs. A scissors 
operator of 27.4 eV is applied to the underlying KS electronic structures to match the experimental absorption onset (Ref. [14]). 
Exciton binding energies, corresponding to the difference between IPA and BSE excitation energies (i.e., Ea — Eaa Eb — Eba 
and Ec — Ec')^ are highlighted in bold. All quantities are expressed in eV. 


crease in presence of CF 3 termination. This functional 
group has an electron-withdrawing character and intro¬ 
duces a sizable dipole moment in the molecule, of the 
order of 5 Debye. This slightly enhances the e-h attrac¬ 
tion, thus strengthening the exciton binding energy of 
the main peaks in Fig. of about 0.2 eV on average. In 
the p-SAM the binding energy of A’ is 0.4 eV larger than 
in the H-az system. On the contrary, in the d-SAM the 
value of Eij decreases by 0.2 eV for the lowest-energy res¬ 
onance A’, compared to its H-terminated counterpart. In 
a similar fashion, the optical absorption onset computed 
for polycyclic aromatic hydrocarbons is red-shifted W 
about 0.3-0.5 eV, in presence of edge-functional groups.l^ 
Also for CFg-az, we observe singlet-triplet splittings (^50 
meV), which are two orders of magnitude smaller than 
the binding energies. The comparison with the experi¬ 
mental data from Ref. [14] indicates good agreement with 
our results. In both theoretical and experimental spectra 
the intense low-energy peak A’, as well as the weaker res¬ 
onances B’ and C’, present the same energy separation 
and relative intensity. Since a* states are not included 
in our BSE calculation, the bump above 405 eV is not 
reproduced in our spectrum. 

IV. SUMMARY AND CONCLUSIONS 

We have investigated N Is x-ray absorption spectra 
of azobenzene-functionalized SAMs, determining the na¬ 
ture of the excitations and discussing the role of many- 
body effects. Our results, obtained from ab initio cal¬ 
culations, in the framework of many-body perturbation 
theory, reveal a clear excitonic character of the main 
peaks in the XAS. Binding energies for core-edge exci- 
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tons, computed from the solution of the Bethe-Salpeter 
equation, decrease from 6 eV in the molecule to 4 eV in 
packed SAMs. This is a many-body effect assigned to an 
interplay between screening and exciton delocalization. 
Based on this finding, we expect exciton coupling be¬ 
tween different molecules to be even more pronounced in 
the optical range, where e-h pairs are typically more delo¬ 
calized. This could give insight into the loss of switching 
capability, as observed for densely packed SAMs.l^^ 

Core-level excitations in these systems are ruled by 
the attractive e-h correlation. The exchange interaction 
plays a negligible role, as testified by singlet-triplet split¬ 
tings, which are two orders of magnitude smaller than 
exciton binding energies. Functionalization with a CF 3 
group does not affect the overall spectral features, but 
only induces a slight increase in the exciton binding ener¬ 
gies of ^ 0.2 eV on average. Good agreement is observed 
with available experimental data. 

In conclusion, our work confirms the predictive power 
of many-body perturbation theory in determining the 
character of the resonances and in disclosing the mi¬ 
croscopic mechanisms ruling the core excitation process. 
This confirms the indispensable role of theory not only 
in interpreting the experimental data, but also in gaining 
further insight into the physics of core-level spectroscopy. 
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